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Abstract 

This paper investigates the performance of open-loop multi-antenna point-to-point links in ad 
hoc networks with slotted ALOHA medium access control (MAC). We consider spatial multiplexing 
transmission with linear maximum ratio combining and zero forcing receivers, as well as orthogonal 
space time block coded transmission. New closed-form expressions are derived for the outage probability, 
throughput and transmission capacity. Our results demonstrate that both the best performing scheme and 
the optimum number of transmit antennas depend on different network parameters, such as the node in- 
tensity and the signal-to-interference-and-noise ratio operating value. We then compare the performance 
to a network consisting of single-antenna devices and an idealized fully centrally coordinated MAC. 
These results show that multi-antenna schemes with a simple decentralized slotted ALOHA MAC can 
outperform even idealized single-antenna networks in various practical scenarios. 
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I. Introduction 

Multiple antennas offer the potential for significant performance improvements in wireless communication systems 
by providing higher data rates and more reliable links. A practical method which can achieve high data rates is to 
employ spatial multiplexing transmission in conjunction with low complexity linear receivers, such as maximum 
ratio combining (MRC) or zero forcing (ZF) receivers. In the context of point-to-point systems operating in the 
absence of interference, the performance of such techniques has now been well studied (see e.g., [1,2]). Multiple 
antennas can also be used to provide more reliable links through spatial diversity techniques. Of the various spatial 
diversity techniques which have been proposed, orthogonal space time block codes (OSTBC) have emerged as one 
of the most important practical approaches, since they offer high diversity gains, whilst requiring only very low 
computational complexity. The performance of OSTBC techniques have also been well studied in the context of 
point-to-point systems, operating in the absence of interference (see e.g., [3,4], and references therein). 

In this paper, we investigate spatial multiplexing techniques and OSTBC in the context of wireless ad hoc 
networks using a slotted ALOHA medium access control (MAC) protocol. The interference model we use includes 
the spatial distribution of nodes, with the nodes distributed as a homogeneous Poisson point process (PPP) on the 
plane. Besides approximating realistic network scenarios, modeling the nodes according to a PPP has the benefit 
of allowing network performance measures to be obtained. This model has been used previously in [5] which 
considered code division multiple access (CDMA) systems, and in [6] which considered single-antenna systems 
with threshold scheduling and power control. It was also used in [7], where network performance measures for a 
single antenna coordinated MAC protocol were derived. 

The PPP model has also been extended to ad hoc networks employing multiple antennas at each transmitting and 
receiving node. In particular, the schemes in [8-11] considered the use of the receive antennas to cancel interference 
from the L closest nodes to the receiver. However, a key requirement of these schemes is that each receiver 
must measure the channels of the L closest interfering nodes, or at least the corresponding short-term covariance 
information, the practical feasibility of which is still not clear. In contrast, the schemes in [12-15] considered the use 
of the receive antennas for either canceling self-interference or increasing the signal power from the corresponding 
transmitter only, and as such, require the receiver to estimate only the transmitter-receiver channel (something which 
can be done with standard channel estimation techniques). For this scenario, [12] considered various spatial diversity 
techniques, and derived transmission capacity and outage probability expressions. In [13, 14], the multiple antennas 
were used for spatial multiplexing, assuming a "closed-loop" scenario where channel knowledge is available at 
the transmitter. In particular, [13] considered multiple-input multiple-output (MIMO) singular value decomposition 
systems, whereas [14] considered MIMO broadcast transmission, for which each transmitting node communicated 
independently with multiple receiver nodes. Note that all of these papers, with the exception of [12, 15], assumed 
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that either the receiver can acquire knowledge about the network interference, or that there is sufficient capabilities 
for each receiver to feedback channel information back to their corresponding transmitter. In practice, it seems 
reasonable to investigate simpler techniques which do not require the receiver to constantly measure the network 
interference, nor require the transmitter to have channel information. We refer to such schemes as point-to-point 
"open-loop" schemes. 

One previous contribution which focuses on such schemes is provided in [15], where OSTBC and spatial 
multiplexing with ZF receivers were considered, and approximations were derived for the frame-error probability. 
In general, however, the fundamental question as to whether to use open-loop spatial multiplexing or diversity-based 
transmission in ad hoc networks, particularly from a capacity-based network performance point of view, is not well 
understood. 

In this paper, we derive new network performance measures for point-to-point open-loop spatial multiplexing with 
MRC and ZF receivers 1 , and OSTBC, using a slotted ALOHA MAC, and the same PPP ad hoc network model as in 
[5, 12]. In all cases, we do not require the receiver to know the network interference, and the transmitter does not have 
any channel state information. (Note that for the ZF case, each receiver cancels interference from its corresponding 
transmitter only.) We derive new closed-form expressions for the outage probability, network throughput, and 
transmission capacity. For spatial multiplexing, these results are exact; whereas, for OSTBC, they are accurate 
approximations, which we show to be significantly more accurate than previous corresponding approximations. 

Our analytical results reveal key insights into the relative throughput performance of spatial multiplexing with 
MRC and ZF receivers, as well as OSTBC, demonstrating that each scheme may outperform the other, depending 
on different network parameters such as the node intensity and the chosen signal-to-interference noise ratio (SINR) 
operating value /3. We also gain insights into the optimal number of transmit antennas to employ for each scheme. 
For certain scenarios, e.g., sufficiently dense networks or networks with high operating SINR values, these insights 
are obtained analytically. For example, we show that single-stream transmission is throughput-optimal in dense 
networks and for networks operating with high (3 requirements, and that transmitting the maximum number of data 
streams is throughput-optimal for systems operating with low requirements. 

We then analyze the transmission capacity of the spatial multiplexing and OSTBC systems. Among other things, 
we prove that for spatial multiplexing, the transmission capacity can scale linearly with the number of antennas, and 

'Note that we consider linear receivers in the spatial multiplexing mode since non-linear receiver structures, such as joint 
decoding, come at the expense of prohibitively higher complexity, especially with increasing number of antennas. We consider 
MRC and ZF receivers, but not the better performing minimum-mean-squared-error (MMSE) receivers, since the analytic 
expressions for MMSE receivers have proved intractable at this stage. However, insights into the MMSE receiver can be directly 
gained from the MRC and ZF results by noting that in both the interference and non-interference scenario, the performance of 
MRC and ZF receivers converges to the performance of MMSE receivers at low and high signal-to-noise ratios respectively. 
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we derive precise conditions which must be met to achieve this scaling, using tools from large-dimensional random 
matrix theory. We also provide concrete design guidelines for selecting the number of transmit streams for achieving 
(and maximizing) this scaling, and contrast our results with those in [10], which derived a similar scaling result 
with receivers employing (network) interference cancelation. In addition, we prove that the transmission capacity 
of OSTBC scales only sub-linearly with the number of antennas, regardless of the particular code employed. 

Finally, we turn from the analysis and comparison of different MIMO schemes to look closer at the benefits of 
multiple antennas with a decentralized MAC, compared with the benefits of a tightly coordinated MAC. In particular, 
we compare the performance of a multi-antenna system employing a simple slotted ALOHA MAC protocol with a 
baseline scheme involving single-antenna devices and a fully coordinated access (CA) MAC which enforces guard 
zones around each receiver. This scheme is similar to that proposed in [7], however we also consider a time division 
multiple access (TDMA) scheme where only one transmitting node around each receiver within the guard zone is 
scheduled to transmit. It is worth noting that this is an idealized protocol, since the overhead involved in achieving 
full coordination is prohibitive in practice for ad hoc networks. To compare with this baseline scheme, we first derive 
new throughput expressions for the CA MAC protocol. These expressions allow us to show the important result 
that the slotted ALOHA approach with multiple antennas can actually perform better than the idealized CA MAC, 
for a wide range of system parameters. This shows that not only does the use of multiple antennas compensate 
for the inherent performance loss caused by using a simple decentralized MAC (i.e., slotted ALOHA), but it can 
actually yield an overall performance improvement. For sparse networks, we show that this is true for most slotted 
ALOHA transmission probabilities, and for dense networks when there is a sufficient number of antennas. 

II. System Model 

We consider an ad hoc network comprising of transmitter-receiver pairs, where each transmitter communicates 
to its corresponding receiver in a point-to-point manner, treating all other transmissions as interference. In addition, 
each transmitter-receiver pair are separated by a distance r tr meters. The transmitting nodes are distributed spatially 
according to a homogeneous PPP of intensity A nodes per meter squared in R 2 . Each transmitting node transmits with 
probability p according to a slotted ALOHA MAC protocol. As such, the effective intensity of actual transmitting 
nodes is pX. 

In this paper, we investigate network performance measures. To obtain such measures, it is sufficient to focus on 
a typical transmitter-receiver pair, with the typical receiver located at the origin. This can be done due to the Palm 
probabilities of a PPP, which states that conditioning on a typical receiver located at the origin does not affect the 
statistics of the rest of the process [16]. In addition, the stationarity property of the PPP indicates that the statistics 
of the signal received at the typical receiver is the same as for every other receiver. Note that the typical transmitter, 
i.e., the transmitter associated with the typical transmitter-receiver pair, is not considered part of the PPP. 
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We consider a network where each node is equipped with TV antennas. Each transmitting node communicates using 
M out of its N antennas, whereas each receiver operates employing all N antennas. As we will show, transmitting 
with less than N antennas may, in fact, lead to an increased network throughput due to the reduced interference in 
the network. The transmitting nodes, with the exception of the typical transmitter, constitute a marked PPP. This 
is denoted by $ = {(De, ~Rt),& € 1N + }, where Di and 2 ~ CN 'n,m (Onxm , Ijvxm) model the location and 
channel matrix respectively of the £th transmitting node with respect to (w.r.t.) the typical receiver. Further, we 
denote the typical transmitter as the Oth transmitting node, and the channel matrix of the typical transmitter-receiver 
pair given by Ho ~ CNn,m (Onxm, Ijvxm)- Each transmitting node is assumed to use the same transmission 
power P, and the transmitted signals are attenuated by a factor l/r a with distance r where a > 2 is the path loss 
exponent 3 . 

We consider the practical scenario where each receiver has perfect knowledge of the channel to its corresponding 
transmitter, but does not know the channel to the other transmitting nodes. Moreover, we assume an open-loop 
scenario where there are no feedback links between the transmitters and receivers, and as such, all transmitters have 
no channel state information (CSI). This is particularly relevant to ad hoc networks, where obtaining CSI may be 
difficult due to their changing nature. 

In this paper, we consider the use of multiple antennas for either open-loop point-to-point spatial multiplexing 
with linear receivers, or OSTBC. There has been little work, to the author's knowledge, done on analyzing and 
comparing these systems in ad hoc networks. We analyze important network performance measures for these spatial 
multiplexing and OSTBC systems, and draw key insights into their relative performance. 

A. Spatial Multiplexing with Linear Receivers 

For spatial multiplexing transmission, we assume that each transmitting node sends M < N independent data 
streams to its corresponding receiver. Focusing on the kth stream, the received N x 1 signal vector at the typical 
receiver can be written as 



m I : m 



yo 



q=l,q^k D e e<l> V 1 *' q=l 



where h Pi9 is the qth column of H p , x PtQ is the symbol sent from the qth transmit antenna of the pth transmitting 
node satisfying E[|a;p i9 | 2 ] = and n .fc <~ CMn,i {0nx1,N I n ) is the complex additive white Gaussian noise 
(AWGN) vector. To obtain an estimate for x _ k , we consider the use of low complexity MRC and ZF linear receivers. 

? d 

The notation X ~ Y means that X is distributed as Y. 

3 Note that there are more accurate, but more complicated, path loss models particularly suited for dense networks; e.g., 
1/(1 + r a ). We have considered such a model, but have found that it leads to more cumbersome expressions, without changing 
the fundamental insights gained based on the simpler model employed in this paper. 



5 



Revised version. Submitted to IEEE Tr. Info Theory. 



For MRC, the data estimate is formed via xo,k = fc yo,fc, where (-)t denotes conjugate transpose, from which 
the SINR can be written as 



7MRC,0,fc 



A^°ll h 0,fe|| 2 



p sr M l h o,fc h ",q I 2 , _p_ q=1 ll h o.fcH 2 , | 

Mr" 2sq=l,q?k ||h ofc || 2 M \D e \°> +1 



where p = ^- is the transmit signal-to-noise ratio (SNR). Note that 



tt" — llhn fell 2 ~ Gamma ( N, — ) , 



Mr: 
P 



(2) 



M 



Mr\ 



h ^'!£ Gamma (m - 1, P 



Mr! 



and 



£ V A Gamma (jtf 



(3) 



where Gamma(v,#) denotes a gamma random variable with shape parameter v and scale parameter 6, with 
probability density function (p.d.f.) 



fx(x) = 



x v ~ x e 



, x>0 



(4) 



I»0" 

For ZF, since each receiver only knows the CSI of its corresponding transmitter, the ZF weight vector is designed 
to cancel interference due to the other data streams (i.e., other than the one being detected) originating from the 
corresponding transmitter only. The data estimate can thus be written as 



*o,fc = gj.feyo.fe 



(5) 



D e e<£ 



where k is the fcth row of (hJh ) Hq. The SINR for ZF follows from (5), and is given by 



7ZF,0,fe — 



Mr? r [(HjH )-i] 



(HtHo)-!Ht(j:^ e<t ^H^tjHofHjHo)- 



(6) 



m[(hJh )- 



+ 1 



where denotes the (k, fc)th element. Note that 

* Gamma [N-M + l, -£—) 



Mr? 



(7) 



and, from [17, Eq. (2.47)] and [18], it can be shown that 

. [(H^Ho)- 1 ^ (H,Hj) H (hJh )- 



M 



(HjHo)-i 



Gamma ( M, ) 
V M/ 



(8) 
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B. OSTBC 

For OSTBC, N s different symbols are transmitted over r time slots using M < N antennas. Different codes 
have been proposed for OSTBC (see e.g., [19-22]), each being characterized by different values of N s , M and r. 
Associated with each code is a code rate, which is defined by R := 

The transmitted OSTBC M x r code matrix is given by 

X e = J2 Kg A 9 + x l q B i) (9) 

q=l 

where xg q is the qt\\ transmitted symbol of the £th transmitting node, and A q and B 9 are M x r matrices, both 
of which are dependent on the particular code employed. The received N x r signal matrix at the typical receiver 
can be written as 



Y = J — H X + V J t-Jt-H/X/ + N (10) 

where N ~ CNn.t (Onxt, N Ipf XT ) is the AWGN matrix. We assume that the channels H and are constant 
during the r time slots used for transmission. To obtain an expression for the data estimate for the fcth symbol, 
it is convenient to introduce the matrix function ^fe(-) : C Nxt -4- C Nxt , which, for a given input matrix V with 
(p, q)th element v p , q , produces the matrix = <rfc(V) with (p, q)th element 



z k,p,q 



v *p,1 


if 




~ v p,q 


if 


v p,q 


if 



(P. l) € ft 



(p,q) G Q\(<P k ux fe ) 
where (•)* denotes conjugate, and q denotes the entire matrix index set given by 

e = {(p,?):pei,... 1 JVn«6i r} , (li) 

with (p C q and x C q. Here, the mapping function <Tfc(-)> and index sets ^ and x,, depend, once again, on the 

ft! /v /v /v 

specific OSTBC code employed. Given these code-specific parameters, the data estimate for the fcth symbol is then 
obtained via the following operation 

x , k = ||M fe ft (Y )||i (12) 

where || • ||i denotes 1-norm (i.e., the sum of all entries), denotes Hadamard product (i.e., elementwise product), 
and M fe is an N x r matrix chosen according to the principle of MRC to satisfy 

||M fc Oft(H Xo)||i = ||H ||^o, fc (13) 

where || • \ \p denotes Frobenius norm. 
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To further illustrate the code-specific parameters A q , B 9 , X^, Qj(-) and M&, in the general OSTBC model 
presented above, let us consider the following concrete example. 
Example: [Alamouti code with N — 2]: In this case, we have 



Ai 



A 2 - 



1 





1 



and B 2 = 




1 

-1 




(14) 



Substituting (14) into (9), the OSTBC code matrix can be written as 

X i,2 X £,l 

Let us focus on decoding x ^ (i.e., k — 1). Then the mapping function ci(-) is given by 



(15) 



ft(V) 



(16) 



"1,1 V l,2 

v 2 .i v* 22 

Substituting (15) and (16) with I = into (13), we find that the matrix Mi which solves the resulting equation is 
given by 

^0,1,1 ^0,1,2 
^0,2,1 ^0,2,2 



Mi = 



(17) 



□ 

Returning to the general case, by noting that ffc(-) and are both linear functions, the data estimate (12) for the 
fcth symbol can be written as 

l H o|lF£0,fc | 



from which the SINR is obtained as 



\ r?Th l|Mfe ft ( H ^) I' 1 + ll Mfe * ( N o) id 

D e e<s> V ' £ ' 



7OSTBC,0,fe = 



Tuk? r II h oIIf 



RM Ed £ 6$ |D*|<*^,£ + 1 



where 



K 



Ex, [||M fc 0<r fc (H £ X £ )||2] 



l H o||| 



(18) 



(19) 



(20) 



is the normalized interference power for the £th transmitting node. 
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Deriving exact closed-form expressions for the SINR distribution is difficult, based on the exact SINR expression 
in (19). To proceed, we focus on deriving an approximation for the SINR, based on assumptions relating to the 
distribution of the JC e _j2 terms, and the independence between different random variables. We explain and justify 
these assumptions in Appendix A, some of which were also used in [12,23]. Our approximation for 7ostbc.o,/« 
derived in Appendix A, is given by 

|Hn 



RMr£ ll"UHF 

TOSTBCAfc = -^—^ 1 5 — " <- 21 ) 

EM Z^D £ e* \D e \" K 'i,'E + 1 



RM 

where 

—tt^r ~ Gamma ( — ) , — - — IIHollp ~ Gamma ( MN, — P —- | (22) 
RM \M'RM) RMr" T UF \ RMrgJ 

are independent random variables, and Nj is the total number of non-zero elements in the columns of the OSTBC 
code matrix containing either ±xe t k or ±x* e k . We will show in Section III that the approximation in (21) is 
very accurate, significantly more so than a previous approximation presented in [12]. 

C. Performance Measures 

In this paper, we consider three main performance measures; namely, outage probability, network throughput, 
and transmission capacity. The outage probability is defined as the probability that the SINR falls below a certain 
threshold (3, i.e. 4 , 

F(f3) = Pr(SINR < /?) . (23) 

The network throughput is defined as the total number of successful transmitted symbols/channel use/unit area. 
This is given by 

T = G>A(1-F(0)) (24) 

where £ is the average number of transmitted symbols per node per channel use, and is given in Table I for spatial 
multiplexing and OSTBC. 

Although throughput is an important performance measure, it may be obtained at the expense of unacceptably 
high outage levels, which is undesirable for some applications due to, for example, significant delays caused by 
data retransmission. This has motivated the introduction of the transmission capacity [5], defined as the maximum 
throughput subject to an outage constraint e, where the maximization is performed over all intensities Xp. The 
transmission capacity is thus given by 

c(e) = CA(e)(l-e) (25) 

4 Note that an alternative outage definition, adopted in [11], would be to consider the probability that the overall sum-rate 
(summed over all data streams) lies below a certain threshold. With this definition, different insights than those presented in this 
paper may be obtained, and this is the subject of ongoing work. 
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TABLE I 

Gamma parameters and £ values for the spatial multiplexing and OSTBC systems. 





SM with MRC receivers 


SM with ZF receivers 


OSTBC 


Desired signal, W ~ Gamma(m, 9) 
Multi-node interference, ~ Gamma(n, Q) 
Self-interference, Y ~ Gamma(u, T) 
Effective multiplexing gain, (" 


m = N, 9 = ^ 
n = M, a = £ 
u = M-l, T= -r/^ 

M 


m = JV-M + 1, ^=m% 
n = M, n = £ 
N/A 
M 


™ = MN >° = Tuk* 
n _ Nj_ n - _e_ 

N/A 



where A(e) is the contention density, defined as the inverse of e = F(/3; Ap) taken w.r.t. Xp. Note that here we have 
made explicit the dependence of the outage probability on Xp. In the following three sections, we will investigate 
each of these three performance measures for spatial multiplexing and OSTBC systems. 

III. Outage Probability and Network Throughput: Exact Analysis 

In this section, we derive new exact closed-form expressions for the outage probability and network throughput 
for the spatial multiplexing and OSTBC systems. To facilitate the derivations, we first note that the received SINR 
for the spatial multiplexing and OSTBC systems can be written in the general form 

W 

7 = Y + Z t€ *\Xi\- a *t + l (26) 

with the generic random variables 

W ~ Gamma (m, 9) , Y ~ Gamma (u, T) ^ m ~ Gamma (n, 0) . (27) 

Here W represents the effective signal power from the desired transmitter, Y represents the effective self-interference 
power from the desired transmitter, and represents the effective interference power from the interfering trans- 
mitting nodes. These are all mutually independent. The particularizations of the shape and scale parameters of W, 
Y, and ^>£o f° r the spatial multiplexing and OSTBC systems are summarized in Table I. 

From Table I, we can directly compare the SINR distributions of each system and make some important 
observations. Focusing on the spatial multiplexing systems, it is evident that the distribution of the effective 
interference <J^o caused by the interfering transmitting nodes is the same for both MRC and ZF receivers. The 
difference lies in the effective signal power W and the self-interference Y. For ZF receivers, the receive d.o.f. 
is split between boosting the signal power and canceling the total self-interference. For MRC receivers, the total 
receive d.o.f. is used to boost the signal power, but no effort is made to cancel self -interference.. 

Based on (26), we present the following general theorem which, after substituting the parameters in Table I, 
yields exact closed-form expressions for the outage probability of the spatial multiplexing and OSTBC schemes. 
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Theorem 1: If the SINR 7 takes the general form (26), then its cumulative distribution function (c.d.f.) is 



F 7 (/3) = 1 



( _ ir -i e -A(4P)^(«) e -§ ^ (m _ I 



T(m) 



E 

£=0 



< m-C-l 



s(m-£,i + i) 



3=0 



1 t 



r=0 



i=0 



j3V 



where 



rj{n) := 



^r(n + f)r(i-f) 

I» 



(28) 



(29) 



and s(n,m) and S(n,m) are Stirling numbers of the first and second kind respectively [24, pp. 824]. Further, 



captures the effects of the self-interference on the outage probability, and is given by 

T(t + u) {(3 1 



e' — Y 1 



(30) 



T(u)T u \0 T, 

Proof: See Appendix B. ■ 
Note that the outage probability and transmission capacity of systems with SINRs of the general form (26) 
were also considered previously in [6] and [12], focusing specifically on the special case Y — (i.e., no self- 
interference). However, in contrast to our results, [6] presented bounds, rather than exact expressions; whereas the 
results in [12], whilst exact, were derived using different methods and were expressed in a more complicated form 
involving summations over subsets. 

Corollary 1: For m = 1 and Y = 0, the SINR c.d.f. (28) becomes 



F 7 (/3) = l- e - A (^)^(») e -§ . 



(31) 



Note that this very simple expression can be used to give the outage probability of spatial multiplexing with ZF 
receivers when M = N. 



A. Throughput of Spatial Multiplexing 

To compute the throughput achieved by spatial multiplexing with MRC and ZF receivers, we substitute the 
relevant parameters from Table I into (28), and substitute the resulting expression into (24). 

Figs. 1 and 2 show the throughputs achieved by spatial multiplexing with MRC and ZF receivers for different 
SINR operating values /3, and different numbers of transmission streams M. In both cases, results are shown for 
N = 4 antennas. We see that for all curves, the throughput increases monotonically with A up to a certain value, 
after which it decreases monotonically. This behavior is intuitive, since increasing A results in a higher number 
of transmissions in the network, however it also yields more interference for a given link. The inherent trade-off 
between these competing factors is clear from the figures. 
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0.04 




0.06 



Intensity, X nodes/m 



Fig. 1. Throughput vs intensity of slotted ALOHA with spatial multiplexing and MRC receivers, and with N = 4, a = 3.1, 
r tr = 2 m, p — 25 dB and p = 1. 



For both receiver structures, we see that in many cases, the throughput is significantly higher when less antennas 
are used for transmission. This is particularly significant if the SINR operating value j3 is high, and when the spatial 
intensity A is large. This behavior can be explained by noting that in these regimes, the additional interference in 
the network caused by each transmitting node employing more antennas outweighs the benefits of an increased 
multiplexing gain. As we may expect, however, for the contrasting scenario where j3 and A are low, it is beneficial 
(in terms of network throughput) to use more transmit antennas. The results in Figs. 1 and 2 demonstrate various 
important tradeoffs which arise in multi-antenna ad hoc networks when using spatial multiplexing transmission. We 
will examine these further in the following section, where we derive simplified expressions for various asymptotic 
regimes. 

1 ) Special Case: Spatial Multiplexing with ZF (M = N): For this special case, as evident from 
(31), the exact outage probability and throughput expressions reduce to very simple forms, and we can gain some 
interesting analytical insights. In particular, taking the derivative of Tzf, we find that with all other parameters 
fixed, the throughput is maximized if A is chosen to satisfy: 



A opt 



T(N) 



(32) 



This point is analogous to the "peaks" identified in Figs. 1 and 2, and gives the optimal tradeoff in terms of network 
interference and spatial multiplexing gain. As we may expect, we see that this optimal spatial intensity decreases 
when N increases (since there is more network interference, whilst the per-link multiplexing gain is unchanged), 
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Fig. 2. Throughput vs intensity of slotted ALOHA with spatial multiplexing and ZF receivers, and with N = 4, a = 3.1, 
r tr — 2 m, p — 25 dB and p = 1. 



or when /3 or r tr increases (since this puts higher reliability requirements on the per-link performance, whilst not 
changing the network interference). Note also that the probability of successful transmission, 1 — F 7zF (/3), obtained 
by using A opt , is e _1 . This is interesting since it coincides precisely with the probability of successful transmission 
for single antenna systems, presented in [6,25]. 

In addition to considering the optimal A, it is also of interest to study the optimal number of antennas. In 
particular, if all other parameters are kept fixed, then we find that the throughput is maximized by choosing N as 

iV opt =max(La;J,l) (33) 

with x the solution to 

1x1(1 + 1)^ + 1) = 2vp\T(l-l)^rl frt 
T(x+l) a p 

This is illustrated in Figs. 3, which plots the optimal number of antennas for different node intensities. We see that 
N opt is decreasing in A, in line with the observations in Figs. 1 and 2. 

B. Throughput of OSTBC 

To compute the throughput achieved by OSTBC, we substitute the relevant parameters from Table I into (28), 
and substitute the resulting expression into (24). Recall that the parameters from Table I give approximations for 
the outage probability and throughput of OSTBC, rather than exact results. 
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Fig. 3. Optimal number of antennas vs. intensity with r tr = 3 m, p — 25 dB, /3 = dB, a = 4 and p = 1. 



Figs. 4, 5, and 6 plot the outage probability vs. SINR threshold for three different OSTBC codes, with N = M. 
Fig. 4 is based on the code 



(35) 



Fig. 5 is based on the code 



and Fig. 6 is based on the code 
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(36) 



(37) 



The analytical curves are seen to accurately approximate the Monte Carlo simulated curves for all SINR thresholds. 
For further comparison, the outage probability approximation given in [12] is also shown. The increased accuracy 
of our approximation is clearly evident. 
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Fig. 4. Outage probability vs SINR of slotted ALOHA with OSTBC using the code in (35), and with a = 3.5, r tr = 5 m, 
A = 0.1 nodes/m 2 , N = 4 and p = 1. 




)• 1 ' 1 ' ' 

0.1 0.2 0.3 0.4 0.5 

SINR Operating Value, P 



Fig. 5. Outage probability vs SINR of slotted ALOHA with OSTBC using the code in (36), and with a = 3.5, r tr = 5 m, 
A = 0.05 nodes/m 2 , N = 3 and p = 1. 



Note that our outage probability approximation becomes an exact result for the important class of cyclic antenna 
diversity codes, which transmit only one symbol per OSTBC codeword, with the symbol being sent out of a different 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



SINR Operating Value, P 



Fig. 6. Outage probability vs SINR of slotted ALOHA with OSTBC using the code in (37), and with a = 3.5, r tr = 5 m, 
A = 0.1 nodes/m 2 , N = 4 and p = 1. 

antenna during each channel use. For example, for the case M = 4, the codeword matrix for this coding scheme 
would have the form 

X, = Xt ,xU • (38) 

This type of code, whilst achieving full spatial diversity order, results in the lowest code rate among all OSTBC 
codes, under the assumption that at least one symbol is transmitted per time slot. However, as we show in Section 
IV, if the network is sufficiently dense, then this type of coding scheme becomes optimal in terms of maximizing 
throughput, due to the minimal network interference it yields compared with other higher rate codes. 

Fig. 7 shows the throughput achieved by OSTBC, based on (24), for different SINR operating values (3, and 
different numbers of antennas used for transmission M. The codes used for M = 2 and M = 3 are given by (15) 
and (36) respectively. Note that both of these codes correspond to maximum-rate codes for the particular antenna 
configuration employed [22]. As with spatial multiplexing, we see that for all curves, the throughput increases 
monotonically with A up to a certain peak point, after which the throughput decreases monotonically. 

Fig. 7 also reveals the interesting fact that the throughput can be significantly higher if less antennas are used for 
transmission. This can be explained by first noting that for M — 2 when using the code in (15), we have ^ = 2, 
and for M — 3 when using the code in (36), we have ^ = | . Hence, whilst increasing M increases the spatial 
diversity order, it also increases (Note that for maximum-rate codes, such as the ones used in Fig. 7, it can 
easily be shown that Nj is always an increasing function of M.) Further, by invoking Lemma 4, we see that the 
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Fig. 7. Throughput vs intensity of slotted ALOHA with OSTBC, and with N = 4, a = 3.1, r tr = 2 m, p = 25 dB and p = 1. 



moments of the approximate normalized interference term ICg in (21) are increasing functions of This suggests 
that the additional interference in the network caused by each transmitting node employing M = 3 antennas using 
the code in (36), compared with M = 2 antennas using the code in (15), outweighs the benefits in terms of an 
increased spatial diversity order. 

Fig. 8 shows the throughput achieved by OSTBC, based on (24), for different SINR operating values (3, and 
different values of TYj. In particular, for TYj = 4, TVj = 8, and TYj = 12, we use the codes in (38), (37), and (35) 
respectively. Note that each of these codes have the same diversity order, but different code rates R. Specifically, 
for the case TYj = 4, R = \; for the case /Vj = 8, R = \; and for the case TYj = 12, R = §. Clearly the code 
rate is an increasing function of Nj. Fig. 8 reveals the interesting fact that in many cases, the throughput is higher 
when a lower code rate is used. This is because the additional interference caused by more simultaneous symbol 
transmissions in the network outweighs the benefits of an increased code rate. 



C. Comparison 

We see in Table I that the signal and interference powers for the spatial multiplexing and OSTBC schemes are 
different, and are dependent on the system parameters M, N, R and Nj. As such, it is not straightforward to 
determine which scheme performs the best, and under which scenario. For example, when M = 2, it can be shown 
that the signal power of OSTBC using the Alamouti code is greater than spatial multiplexing with ZF receivers, 
while the interference power of both OSTBC and spatial multiplexing schemes are the same. However, for a fixed 
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Fig. 8. Throughput vs intensity of slotted ALOHA with OSTBC, and with N = 4, a = 3.1, r tr = 3 m, p = 20 dB and p = 1. 

/3, the per-link data rate for spatial multiplexing is twice that of OSTBC, due to the M — 2 transmitted streams 
per node for spatial multiplexing. Thus although the SINR of OSTBC is greater than spatial multiplexing with ZF 
receivers for M = 2, the overall network throughput may not be greater. 

The best performing receiver used for spatial multiplexing is also dependent on the particular scenario. For 
example, it is well known that for a single user MIMO system with no interfering nodes, spatial multiplexing 
with ZF receivers performs better than spatial multiplexing with MRC receivers at high SNR, hence ZF receivers 
is expected to perform better than MRC receivers in sparse networks. However, for MRC receivers, as the node 
density increases, the impact of the interference from other transmitting nodes becomes more dominant than the 
self-interference. By noting that the distribution of the interference power from these interfering transmitting nodes 
for both MRC and ZF receivers are equal, the distribution of the total interference for both MRC and ZF receivers 
thus converge with increasing node density. The key factor which distinguishes between these two schemes in dense 
networks is the signal power, which is greater for MRC receivers. This suggests that the best performing receiver 
used for spatial multiplexing is dependent on the particular node density. 

The preceding discussion suggests that the best performing scheme is dependent on the system and network 
parameters. This can be seen in Fig. 9, which plots the network throughput vs. node intensity A, of spatial 
multiplexing with MRC and ZF receivers, and OSTBC, based on (24), for different numbers of antennas used for 
transmission, M. We observe that spatial multiplexing with ZF performs the best in sparse network configurations 
(i.e., small A), while OSTBC performs the best in dense networks (i.e., high A). Further, spatial multiplexing with 



18 



Revised version. Submitted to IEEE Tr. Info Theory. 




Fig. 9. Throughput vs intensity of slotted ALOHA with spatial multiplexing and OSTBC, and with a = 3.3, r tr = 2 m, 
p = 25 dB, /3 = 2 dB and p = 1. 

MRC receivers is seen to perform better than ZF receivers at high A, which agrees with previous discussion. We 
will also prove this analytically in Section IV. 

The best performing scheme can also be shown to be dependent on the SINR (3 operating value, as shown in 
Fig. 10, which plots the network throughput vs. SINR (3 for all three schemes, again for different M. We see that 
spatial multiplexing with MRC receivers performs the best for low /?, while OSTBC performs the best for high /3. 

As can be seen from Figs. 1-10, the impact of the number of antennas used for transmission, and the relative 
throughput of the spatial multiplexing and OSTBC systems are dependent on different network parameters. To 
explore this further, it is convenient to analyze the throughput of spatial multiplexing and OSTBC in asymptotic 
regimes. This is considered in the following section. 

IV. Network Throughput: Asymptotic Analysis 

In this section, we analyze the performance of spatial multiplexing and OSTBC in dense networks, and for low 
and high f3 operating values. We note that in sparse networks, the performance of spatial multiplexing and OSTBC 
approaches the performance of single user MIMO systems, for which the performance is well known (see e.g., [1- 
4]). Thus, although we investigate dense network scenarios, we do not consider the opposite case of sparse networks 
in this paper. Note also, that for M = 1, corresponding to single-input multiple-output (SIMO) transmission, the 
performance of all schemes are the same. 
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Fig. 10. Throughput vs SINR operating value /3 of slotted ALOHA with spatial multiplexing and OSTBC, and with A = 4, 

a = 2.5, r u = 3 m, p = 20 dB, A = 0.01 nodes/m 2 and p = 1. 



A. Dense Networks (Large X) 

Under dense network conditions, we have the following corollary. 
Corollary 2: As A -> oo, the SINR c.d.f. in (28) behaves as 5 

F 7 , A ^co(/3)+of e - A (^)^A m - 



F 7 G8) 



(39) 



where 



F 7 ,A^co(/3) = 1 



277(71) //3fi 



and 



ar(m) V 



e e e 



(40) 



k T/3 + fl y 

To give an indication as to when a network is "sufficiently dense" such that the expansion (40) is accurate, Table 
II tabulates the quantity 1/A m j n , with A m ; n representing the minimum node density required such that (40) is within 
at least 85% of the true non-asymptotic value. Intuitively, the quantity 1/A m ; n gives a measure of the maximum 
allowable separation between adjacent transmitting nodes (on average), in order for the asymptotic expansion (40) to 
serve as a good approximation. Note that a scenario with M = 1 data streams is chosen, because as we will discuss 

5 The notation f(x) = o(g(x)) as x — > 00 means that for every M > 0, there exists a constant xq such that |/(a;)| < M\g(x)\ 
for all x > xq. 
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TABLE II 

Average maximum area around each transmitting node, such that 



3 (/3)-F 7 (« 



1-F 7 03) 

SINR THRESHOLDS j3 AND TRANSMIT SNRS p WITH M = 1, N = 3, r t r = 5 M AND Q = 2.1. 



> 0.15, FOR VARIOUS 



p dB \/3 dB 


1 


2 


5 


10 


20 


2 


6.500 


6.486 


6.383 


6.114 


5.506 


5 


12.447 


12.519 


12.520 


12.133 


10.981 


10 


33.478 


34.674 


36.900 


37.481 


34.638 


20 


111.483 


130.208 


193.050 


286.532 


335.569 



later, single-stream transmission is throughput-optimal in dense networks. Moreover, a transmit-receive distance of 
5 m is chosen, which is practically relevant (e.g., for wireless local area networks). 

The results in the table show that the large-A expansion (40) serves as an accurate performance measure for 
practical network configurations; in some cases, applying even when the networks are relatively sparse. This is 
particularly true for moderate to large SNRs. For example, at 10 dB SNR, a transmitter spacing of roughly 35 m 2 
on average is sufficient, which is relatively large compared with the transmit-receive distance of 5 m. As the SNR 
is reduced, e.g., to 2 dB, a closer transmitter spacing of roughly 6 m 2 is needed. This behavior is intuitive, since 
dense networking conditions are representative of interference-limited scenarios, in which case the interference in 
the network is much more significant compared with the noise. If the SNR is reduced, then the noise has greater 
relative effect, and there must be more interference (i.e., a greater density A) in order for the interference-limited 
behavior to be apparent. 

From (24) and (39), and recalling that Q./6 = r" T for each scheme (c.f. Table I), it follows that for large A the 
throughput becomes 

Tae -A^r? r ,W A m _ (41) 

From this, we can obtain some useful insights into the network performance and optimization: 

• The throughput decays exponentially in A. This loss of throughput indicates that in dense networks, the 
negative effects of interference will dominate any positive throughput gains obtained by an increase in the 
number of communication links. 

• Since the exponential in (41) dominates for large A, the component 77(71) is a critical factor which determines 
performance. Moreover, since n is proportional to the effective interference power caused by each interfering 
transmitter, and 77(71) increases with n, for all three transmission schemes it is best to choose n as small as 
possible. Thus, recalling the parameters in Table I, we have the following design criteria for dense networks: 
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- For spatial multiplexing with either MRC or ZF receivers, it is optimal to use only a single transmit 
stream (i.e., M = 1). 

- For OSTBC, it is optimal to use a cyclic antenna diversity coding scheme. This code, illustrated in 
(38) for M = 4, minimizes n by maximizing the coding parameter Nj (i.e., Nj — M). To determine 
the optimal M for cyclic antenna diversity codes, we consider the leading order factor A m . As m is 
proportional to the spatial diversity order, we note that the optimal choice is M = N. 

• For systems with 6 M > 1, OSTBC codes with Nj < M 2 will yield a higher throughput than both spatial 
multiplexing schemes; whereas, if Nj > M 2 , the throughput will be worse than both. This can be seen by 
again noting that the exponential in (41) dominates for large A. As Nj = M for cyclic antenna diversity code, 
we see that these codes always perform better than spatial multiplexing. For the remaining case Nj = M 2 , 
which only occurs for the Alamouti code and therefore M = 2, the throughputs of all three schemes have 
the same exponential decay in (41); however, OSTBC has the largest leading order factor A™\ and therefore 
achieves the highest throughput. Finally, for all M > 1, spatial multiplexing with MRC achieves a higher 
throughput than ZF. 

In general, these results reveal that OSTBC with a M = N cyclic antenna diversity code is the optimal scheme in 
dense networks. 



B. Networks with High SINR Operating Values (Large (3) 

For networks with high SINR operating values, we have the following corollary. 
Corollary 3: As (3 -> oo, the SINR c.d.f. (28) behaves as 

{J3) + o( E Y e-h-^^^T- 



F 7 (/3) 



■ 7,/3— too 



where 



? 7 ,^ 00 (/3) = l- r(m)6)m _ 1 E y 



(42) 



(43) 



To check the accuracy of this expansion, Table III tabulates for various scenarios the minimum SINR operating 
value j3 such that (43) is within 85% of the true non-asymptotic value. These results demonstrate that the large-/? 
expansion (43) serves as an accurate performance measure for practical network configurations; in some cases, 
applying even when j3 is low. This is particularly true for low to moderate SNRs. For example, at 2 dB SNR and 
A chosen such that the average transmitter spacing is 20 m 2 , a SINR operating value of —12.832 dB or above is 
sufficient — something which is expected to be true for most wireless applications. As the SNR is increased, e.g., 

6 Although M > 1 is sub-optimal for dense networks, for other networking scenarios this configuration will become important. 
Thus, it is still of interest to study the throughput of such configurations under dense network conditions. 
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TABLE III 

F 7l3 - ) .oo(«-F 7 ( / 3) 



Minimum SINR f3 (dB) operating value such that 



l-F T (/3) 

SNRs p WITH M = 1, N = 3, r t r = 5 M AND a = 4. 



> 0.15, FOR VARIOUS 1/A, AND TRANSMIT 



p dB \±m 2 


100 


50 


20 


12.5 


10 


2 


-14.547 


-14.078 


-12.832 


-11.543 


-10.701 


5 


-11.421 


-10.804 


-8.925 


-7.211 


-6.142 


10 


-5.933 


-4.827 


-1.726 


0.885 


2.366 


15 


-0.078 


1.833 


6.569 


9.948 


11.699 



to 10 dB, a higher SINR operating value of —1.726 dB or above is needed. This matches with intuition, since, if 
the network is relatively sparse (as in the example above), then for the outage probability to remain constant as the 
SINR is increased, the SNR must increase accordingly. 

From (24) and (42), and recalling that ft/9 = r" T for each scheme (c.f. Table I), it follows that for large /? the 
throughput becomes 

1 



T oc 



-e o e 



f p -Art r /3 a 7j( n ) am-l-u 



(44) 



r(m)0 m - 1 T> 

From this, we can obtain some useful insights into the network performance and optimization: 

• Recalling that a > 2, we see that the exponential in (44) dominates for large j3. This is intuitive, by 
recalling the direct correspondence between the SINR operating value (3 and the outage probability. 

• The component 9, which represents the transmit SNR, is a critical factor, and this quantity should be maximized 
in order to maximize the throughput. When 9 is independent of M, as is the case for cyclic antenna diversity 
codes, the fj, m - x ~ u polynomial term should be maximized in order to maximize the throughput. Recalling 
the parameters from Table I, we arrive at the following design criteria for networks with high SINR operating 
values: 

- The optimal transmission scheme is the same as for dense networks. That is, for spatial multiplexing 
with either MRC or ZF receivers it is optimal to use only a single transmit stream, whereas for OSTBC 
it is optimal to use the M = N cyclic antenna diversity coding scheme. 

Note that if the SNR is also sufficiently high, then the exponential e _Ar tr' 9 " ifa) m ay dominate e~ ^ ; however, 

it is easy to see that the same optimality criteria still applies. 

• For systems with M > 1, OSTBC codes yield a higher throughput than each of the spatial multiplexing 
schemes. Moreover, since the throughput of both spatial multiplexing schemes have the same exponential 
decay and also the same polynomial factor ^J™- 1 -" i n (44), their relative performance is determined by the 
constant factor r ( m )8™- 1 T" • Thus, by substituting the relevant parameters from Table I, we can show that ZF 
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Maximum SINR f} (dB) operating value such that 



TABLE IV 



I-F7O 3 ) 



> 0.15, FOR VARIOUS 1/A AND TRANSMIT 



SNRs p. Spatial multiplexing with MRC receivers is considered with M = 4,JV = 4,r t r = 5M and a = 3. 



p dB \\m 2 


100 


50 


20 


12.5 


10 


10 


-22.055 


-22.660 


-24.559 


-27.033 


-30.362 


15 


-18.164 


-19.462 


-22.518 


-25.850 


-29.788 


20 


-15.476 


-17.467 


-21.487 


-25.346 


-29.586 


25 


-14.0351 


-16.501 


-21.068 


-25.171 


-29.547 


30 


-13.435 


-16.128 


-20.921 


-25.100 


-29.508 



will deliver a higher throughput than MRC if the following condition is met: 

i 

-1 

(45) 



otherwise MRC will perform better. 
In general, these results indicate that, as for dense networks, OSTBC with a M = N cyclic antenna diversity code 
is the preferable scheme for high SINR operating values. 

C. Networks with Low SINR Operating Values (Small (3) 

For networks with low SINR operating values, we have the following corollary. 
Corollary 4: As /3 ->• 0+, the SINR c.d.f. (28) behaves as 7 

F 7 (/?) = F 7i/? ^ +(/?) + C>(/3) (46) 

where 



r (m-|) V (n) (0Q 

r(m) r(i-f) v 



To give an indication of when the expansion (47) is accurate, Table IV tabulates the maximum SINR operating 
values such that (47) is within 85% of the true non-asymptotic value. From the table, we see that quite low SINR 
operating values are required, and these may or may not be practical. Nevertheless, as we discuss shortly, the 
insights we will obtain from (47) also turn out to be valid for much higher SINR operating values than those 
tabulated in Table IV. 

7 The notation f(x) = 0(g(x)) as x — > + means there exists positive numbers 6 and M such that \f(x)\ < M\g(x)\ for 
\x\ < S. 
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From (24) and (46), and recalling that Q./6 = r" T for each scheme (c.f. Table I), it follows that for small (3 the 
throughput becomes 



^^^- "W^H) ^'] • (48) 

From the preceding equations, we can obtain the following useful insights: 

• If (3 is "very" small (i.e., such that the /3-dependent term in (48) is negligible), then the throughput approaches 
(Xp, which trivially gives the following design criteria: 

- For spatial multiplexing with either MRC or ZF receivers, it is optimal to use the maximum number of 
transmit streams. 

- For OSTBC, it is optimal to use maximum-rate codes. This follows by noting that the maximum code 
rate is a decreasing function of the number of antennas used for transmission when M > 2 [22], and 
thus increasing the number of antennas for transmission leads to a lower throughput. The codes with the 
maximum-rate corresponds to either SIMO transmission with M = 1 or the Alamouti code with M = 2. 
From (46), the Alamouti code performs better than SIMO if 

r (*-i) r w / 2 \ 

T{N)T(2N-l) > V + a) ■ ^ 
We observe that this occurs for low path loss exponents, i.e., as a — » 2. 

Intuitively, the outage probability is very close to zero, and therefore it makes sense to send as much 
information as possible during each channel use. For this scenario, spatial multiplexing thus performs better 
than OSTBC. 

• More generally, for all (3 values for which (48) is accurate (i.e., the /3-dependent term in (48) is not necessarily 
negligible), we can derive the following conditions: 

- For spatial multiplexing with MRC and ZF receivers, the optimal number of data streams is given by 

M opt = min (max ( [x\ , 1) , N) (50) 
where, for MRC, x is the solution to 

_J_ = £(tl)£(^i±|) ( l + ^L\ ( 5i) 

Xvpripi T{N) T(x-1) { ^ a(x - 1) J ' K > 

whilst, for ZF, it is the solution to 

i _ r - x + 1 - g) r (.x - 1 + f ) / 2x 2(x-i) 

T(N-x + l) T(x-1) \ a(x - 1) a(N - x + 1) 
Note that the right hand sides of both (51) and (52) are increasing in x, hence M opt is unique. 

- As N ->• oo, 

^ Q° pt (53) 
N 
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where, for MRC, 

n opt _ 

~ /3 VA^pr t 4 > + 2) 



' ( ; ;. I (54) 



whist, for ZF, g opt is the solution to 



f* = IS^^Z ^ . (55) 



1 



These equations confirm the intuition that, for both receivers, as the node density is increased, less transmit 
antennas should be used. This is because adding more nodes to the network increases the aggregate interference, 
and thus, it is better for each node to transmit will less data streams in order to "balance" the overall network 
interference. This behavior was also observed for SINR operating values as high as 3 dB in Figs. 1 and 2. 

• The results above also reveal that less transmit antennas should be used if the SINR operating value /3 
increases, a phenomenon which is confirmed in Fig. 10. Moreover, the optimal number of transmit antennas 
is at least as many for MRC as for ZF. For the example shown in Fig. 10, for MRC, all transmit antennas 
should be used when (3 is below —3 dB; whereas for ZF, f3 must be below —10 dB. 

• By noting that R(M) < 1, it is clear that spatial multiplexing performs better than OSTBC. Moreover, for 
spatial multiplexing, since is decreasing in m, for M > 2, MRC achieves a higher throughput than 
ZF. 

In general, these results indicate that spatial multiplexing with MRC, with all transmit antennas active, is the most 
favorable scheme for low SINR operating values. 

V. Transmission Capacity 

In this section, we turn to the analysis of transmission capacity. In general, we find that for both spatial 
multiplexing and OSTBC, an exact analysis of the transmission capacity is intractable due to the complexity 
involved with inverting the exact outage probability expressions. One exception is the case of spatial multiplexing 
with ZF receivers with M — N, for which an exact expression for the transmission capacity is obtained from (31) 
and (25) as 

fs N(l-e) (. ( 1 \ f3r«N\ 
czf (e) = — — — log - . (56) 

For all other scenarios, we focus on studying the transmission capacity for small outage levels, which is representative 
of practical systems. Under these conditions, based on (26), we present the following key theorem which, after 
substituting the parameters in Table I, yields closed-form expressions for the transmission capacity of the spatial 
multiplexing and OSTBC schemes. 
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Theorem 2: If the SINR takes the general form in (26), then the transmission capacity as e — F su (/3) — > can 
be written as 

( } n^nir („+l)r(m- f ) Ey [e"^^ (l - m, 1 + | - m, ^)] ^ J 

(57) 

where 1F1 (■,-,■) is the Kummer confluent hypergeometric function, and the notation (•)+ implies (a)+ = max(a, 0). 



Also, the expectation in (57) can be written as 

1\ / Ft\ e / 9\ JL /A 

e- — Y 1 



§_ m-l 

fSY _ 



«"'-r&S("7 I )(l) r H-!)SW" 

7SU 



(58) 



and F bu (/3) represents the outage probability with no multi-node interference, computed as 



fe = T=0 v 



(59) 



The remaining expectations are given in closed-form in (30). 

Proof: Follows by taking a Taylor expansion of (28) around A = 0, and then finding the inverse of the resulting 
expression w.r.t. \p. ■ 

Note that the factor F su represents the outage probability of a single user MIMO system, for which outages 
are caused by self-interference and AWGN. Due to the additional multi-node interference in ad hoc networks, any 
specified outage constraint e which falls below F su (/3) can never be met, and therefore the transmission capacity in 
such cases is zero. This phenomenon is illustrated in Fig. 11, where the outage probability is plotted versus intensity 
for spatial multiplexing with MRC receivers, for different antenna configurations. For the results shown, an outage 
probability of less than 0.2 can never be achieved when M > 2, due to the effects of AWGN and self-interference, 
which ensures that e > when A = 0. Note however, that for OSTBC and spatial multiplexing with ZF receivers, 
there is no self-interference, and F su (/3) in this case accounts for outages due to AWGN only. 

To compute the transmission capacity achieved by spatial multiplexing and OSTBC, we substitute the relevant 
parameters from Table I into (57). The accuracy of our transmission capacity expression is confirmed in Fig. 12, 
which plots the transmission capacity vs. outage probability e for spatial multiplexing with ZF receivers. We see 
that for outage probabilities as high as £ = 0.1, our expression is accurate. Although not shown, similar accuracy 
has been observed for the other schemes also. 

For large numbers of antennas, we present the following corollary for the transmission capacity. 

Corollary 5: As M, N — > oo, the transmission capacity as e — F su (/3) -4- satisfies 

Proof: See Appendix C. 
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Fig. 11. Outage probability vs. intensity of slotted ALOHA with spatial multiplexing using MRC receivers, and with N = 4, 
^tr = 3 m, ol = 4.23, j3 = 3 dB, p = 30 dB and p = 1. 
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Fig. 12. Transmission capacity vs. outage probability of slotted ALOHA using ZF receivers, and with r tr = 3 m, j3 — 3 dB, 
a = 4, N = 4 and p = 1. 



(60) implies that if e > F su (/3) and — — ^ n - converges to a constant for large M,N, then the transmission 
capacity will scale linearly with (. In the following, we will use these results to investigate how the transmission 
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capacity scales with the number of antennas, for the spatial multiplexing and OSTBC systems which we consider. 
Under some conditions, we will show that linear scaling is indeed possible. 

A. Transmission Capacity of Spatial Multiplexing 

Here we use Corollary 5 to gain insights into the scaling behavior of the transmission capacity for spatial 
multiplexing with MRC and ZF receivers, as given by the following corollaries. 

Corollary 6: As A — > oo with M = kN where < k < 1, the transmission capacity of spatial multiplexing 
with MRC receivers as e — H {fi — /3mrc) ~> behaves as 

^ - « + l) " ((/Wc ^tf (£ ' g( Vi MRc))+ + O (((e -H((>- M) 2 ) + ) <«> 
where 

/3mrc = — 7-5 r- (62) 

and TJ(-) is the Heaviside step function. 

Proof: Follows by substituting the relevant parameters for MRC receivers in Table I into (60), and noting that 
[26, Eq. (8)] t^jt converges in probability to , i — c- as A — > 00 with M = kJV. ■ 
Corollary 7: As A — > 00 with M = kJV where < k < 1, the transmission capacity of spatial multiplexing 
with ZF receivers as e — H ((3 — /3zf) — > behaves as 

^ - (?) * ((A- - « + ) ! ( "^f Fr + ° ((<« - " iP - .M> 2 ) + ) m 

where 

^zf = 4 f 1 - A ■ (64) 



Proof: Follows by substituting the relevant parameters for ZF receivers in Table I into (60), and noting that 
[26, Theorem 7.2] W converges in probability to ^ — l) as A — > 00 with M = kJV. ■ 
At low outage probabilities e, these results imply that for MRC 8 , 

6 (A) , /?MRC > /3 

CMRc(e) = < 

0, otherwise 

while for ZF, 

czF(e) = < 



(65) 



9(A), p ZF >f3 
0, otherwise 



(66) 



8 The notation g(N) = 0(/(A)) means that (pif(N) < g(N) < fcf{N) for N -> 00, where 0i and <^ 2 are constants 
independent of A. 
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Moreover, from (61) and (63), assuming that /?mrc > P an d Pzf > P, the transmission capacity can be approximated 
for large M and N as 

cmrcW^M 1 -* (n-MP-^*-Y —L-t, (67) 
V P J nrlP* 

c zfW ,m-s U-m-MZS.) 1 '. 

V P ) Trrl(3« 

From the preceding equations, we observe the following: 

• Under the assumption that the operating SINR p is sufficiently small such that /3mrc > P an d Pzf > P, both 
receivers achieve a transmission capacity scaling which is linear in N. 

• ZF receivers are preferable over MRC receivers for sufficiently high SINR operating values, i.e., (3 > 1, and 
vice-versa for small SINR operating values. 

• Spatial multiplexing with MRC receivers achieves the same performance as a spatial multiplexing system 
with ZF receivers, where the ZF receivers utilize N — M(3 + 1 degrees of freedom (d.o.f.) to boost the signal 
power and M(3 — 1 d.o.f. to cancel the self-interference from Mj3 — 1 self-interfering streams. 

This preceding discussion only tells part of the story, and to gain a clearer picture we must also investigate 
the conditions on (3 in (65) and (66) which ultimately dictate when linear scaling is achievable for each receiver. 
Comparing these conditions, we find that 

Pmrc <Pzf 4- > T~ — • < 68 ) 

This result is interesting, since it gives insight into which receiver, MRC or ZF, can achieve linear scaling "more 
easily", i.e., for a larger range of operating SINRs (3, in terms of the average received SNR pr^ a and the antenna 
parameter k. In particular, (68) will be satisfied when either the average received SNR is sufficiently high or n is 
sufficiently small (i.e., M <C N), in which case, ZF will achieve linear scaling over a wider range of operating SINRs 
than MRC. This is intuitive, since by increasing the average received SNR, the SINR of ZF grows proportionately, 
while the SINR of MRC approaches a constant. Moreover, by decreasing k, the signal power of ZF will converge 
to the signal power of MRC, while the self-interference power of MRC approaches a constant. On the other hand, 
if the average received SNR is sufficiently low or k is sufficiently high (i.e., M rts N), then MRC will achieve 
linear scaling for a wider range of j3 values than ZF. 

In general, from a transmission capacity perspective, these results imply that for networks operating with low 
SNR it is highly desirable to employ MRC receivers in favor or ZF. If the SNR is not low, however, the choice 
is less clear, and depends largely on the specific network operating conditions. In particular, if the SNR is large 
enough that (68) is satisfied, then ZF will have the advantage of achieving linear scaling more easily than MRC. 
Regardless of this point, however, if f$ is chosen sufficiently small, then the condition (68) becomes irrelevant, both 
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schemes will achieve linear scaling. Thus, the specific operating SINR f3, as well as the SNR, play a critical role 
in deciding which spatial multiplexing receiver is most advantageous. 

From (65) and (66), we know that if M = kN, then linear scaling of the transmission capacity is achieved for 
both receivers. A natural question which arises, however, is whether one can do even better than linear scaling, if 
M is not confined to vary linearly with N. This is answered in the following proposition: 

Proposition 1: For spatial multiplexing systems with either MRC or ZF receivers, linear scaling of the transmis- 
sion capacity, i.e., c(e) = Q(N), achieved when M — kN with < k < 1, is the best possible scaling. Moreover, 
if M varies sub-linearly with N, then the transmission capacity scales sub-linearly also. 

Proof: See Appendix D. ■ 

Note that this result is certainly not obvious, since in general, the transmission capacity is typically maximized 
by selecting n < 1; i.e., by choosing only a subset of antennas for transmission. Thus, it is not immediately clear 
whether the optimal M should grow in proportion to N, or at some slower (i.e., sub-linear) rate. 

We now consider the optimal antenna configuration which maximizes the transmission capacity scaling, given in 
the following proposition. 

Proposition 2: The antenna configuration which optimizes the transmission capacity is 

M = max (1, min ([iVK*J , N)) (69) 

where 



K = < 



(1 - ^) . * , for MRC 

(l-?-)—^, forZF 

V a.) 1 | r tr£ ' 



(70) 



From this optimal antenna configuration, we observe the following: 

• The optimal number of data streams is higher for MRC than ZF for sufficiently small SINR operating values 
f3, i.e., j3 < 1, and vice-versa for sufficiently high SINR operating values. 

• As a — > 2, k* — > 0, and the optimal number of data streams decreases. This can be explained by noting 
that low path loss exponents correspond to scenarios where all multi-node interference, including far-away 
interference, is significant. In this scenario, the number of data streams should be decreased to reduce the 
impact of multi-node interference. 

• As a goes large, the impact of multi-node interference becomes negligible, and the optimal number of data 
streams is dependent on the specific network operating conditions. For example, when MRC receivers are 
employed, transmitting the maximum number of data streams is optimal for sufficiently low f3, and vice- 
versa for high ji. For ZF receivers, transmitting the maximum number of data streams is optimal for high 
average received SNR pr^ r a . This can be observed in Fig. 13, which shows the transmission capacity of 



31 



Revised version. Submitted to IEEE Tr. Info Theory. 



x 10 



-6 



6 




,| 1 1 1 1 1 1 L_ 

2 2.5 3 3.5 4 4.5 5 5.5 



6 



Path Loss Exponent, a 



Fig. 13. Transmission capacity vs. path loss exponent of slotted ALOHA with spatial multiplexing using ZF receivers in the 
interference-limited regime, and with r tr = 3 m, /3 = 3 dB, e = 0.0001 and p = 1. 

spatial multiplexing with ZF for different M at high transmit SNR p. The "Analytical" curves in the figure 
are based on (57). 

Proposition 2 also provides a useful design criteria for choosing how many data streams should be transmitted 
for a given path loss exponent. For example, for ZF, when a = 4 and at high transmit SNR p, we see in (69) 
that the optimal number of transmitted data streams should be half the number of receive antennas. The optimal 
number of transmit antennas was also considered in [8], which considered the use of MMSE receivers to cancel the 
interference from the closest transmitting nodes and a single-link performance measure. In contrast, we consider 
the use of ZF receivers to cancel interference from the corresponding transmitter only, and a network performance 
measure. 

1 ) Comparison with Multi-Node Interference Cancelation (Single-Input Multiple-Output with 
Partial Zero-Forcing): Linear scaling was also previously shown to occur in [10] for single-stream transmission, 
i.e., M = 1, by employing a partial zero forcing (PZF) scheme, where the receive antennas were used to 
simultaneously cancel the interference from the fcpzp closest interfering nodes and boost the signal power from the 
corresponding transmitter. Specifically, the transmission capacity of the PZF scheme at high transmit SNR p was 
shown to scale as 




(71) 
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Linear scaling was then achieved by setting fcpzF = kN where < k < 1. At high transmit SNR p, comparing 
(67) and (71), we see that setting fcpzF = M results in the PZF scheme achieving the same scaling as our spatial 
multiplexing with ZF receivers scheme 9 . 

To understand the reason behind this same linear scaling result, we explore some similarities between our spatial 
multiplexing with ZF receivers scheme and the PZF scheme in [10]. For both schemes, it can be shown that the 
signal component contributes to the transmission capacity scaling by a factor of N%, and is achieved by allocating 
N(l — k) d.o.f. to boost the signal power, where < k < 1. The key difference, however, lies with how the 
remaining Nk d.o.f. are used. For our spatial multiplexing with ZF receivers scheme, the Nk d.o.f. are used to 
cancel interference from the corresponding transmitting node, while for the PZF scheme, they are used to cancel 
interference from the fcpzF = Nk closest interferers. However, it turns out that the resulting interference component 
of both schemes contributes to the transmission capacity through the same scaling factor, which is N 1 ^^ . To see 
why, it is convenient to express the transmission capacity (25) in the alternative form 

c ZF (e) =A*(e)(l-e) (72) 

where A*(e) is the inverse of e = F ZF taken w.r.t. A. This alternate form can be obtained by noting that 

e = F(/?;A/M) => A = MF" 1 (e), (73) 

which has the same form as (25) for small e. (72) indicates that from a transmission capacity perspective, increasing 
the number of data streams M has the same effect as decreasing the density of transmitting nodes by a factor of 
M. 

From the modified transmission capacity definition (72), we can now explain why for both schemes the interfer- 
ence component contributes to the transmission capacity by the same scaling factor. First, for spatial multiplexing 
with ZF receivers, the average distance from the receiver to the closest interfering node, where the interfering nodes 
are distributed with density A/M, can be shown to scale as \[M [27, Eq. (9)] whereas for PZF (after cancelation), 
it scales as v^pzf [10]. Thus, with M = fcpzF, these coincide. Second, for both schemes, the average number of 
interfering data streams per unit area is the same, given by A. For spatial multiplexing with ZF, this can be seen 
by noting that M data streams are being transmitted by A/M nodes per unit area. 

Summing up, the key point of the discussion above is that the same linear transmission capacity scaling can be 
achieved whether we employ MIMO with ZF receivers, or SIMO with PZF. Clearly, the choice as to which scheme 
to employ in practice will depend on the specific network design. When there is the capability for employing 
multiple antennas on all nodes, MIMO with ZF will be more attractive, since it has lower complexity and less 
stringent requirements on the level of channel knowledge at the receivers. On the other hand, if single-antenna 

9 Note that similar comparisons can be made with our spatial multiplexing with MRC receivers scheme. 
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transmitters are required (e.g., due to size limitations, relevant for sensor networks), then SIMO with PZF will be 
appropriate. A key implementation issue for such networks, however, is how to accurately estimate the multi-node 
interference which is needed to perform interference cancelation at each receiver. Whilst very well motivated in 
theory, some important questions still remain as to the practical feasibility of this approach. 

B. Transmission Capacity of OSTBC 

We now consider the transmission capacity of OSTBC at low outage probabilities e. In this case, analogous to 
(67), we consider the case where N grows large, and also assume that f3 and e are selected such that e > F su (/?). 
Moreover, we consider high transmit SNRs p, and thus, substituting the relevant parameters from Table I into (57), 
we get 

CosTBc(e) - (g(M)Ni) (74) 

where 

r / Nj(M) \ 

g(M) = R(M)Mi V ^ / (75) 

r ( Nj(M) 2 \ 
L \ M a J 

where we have made explicit the dependence of the code rate R and Ni on M. (Note that this result applies for 
arbitrary values of M.) In general, there is no simple relationship characterizing the dependence of R and Nj on 
M. Thus, to proceed, here we consider the extreme scenarios pertaining to minimum-rate and maximum-rate codes. 

For minimum-rate codes, i.e., cyclic diversity systems, R(M) = 1/M and Nj(M) = M. Thus, it follows that 
Costbc( £ ) = @ (^Mi~ 1 Ni^J, which is obviously maximized if M = 1, yielding c stbc( £ ) = © (^N*^. This 
special case simply corresponds to a standard SIMO system. 

For maximum-rate codes, e.g., the Alamouti code, it is known from [22] that for M even 

and for M odd 

Thus, using these bounds and the inequality < T \^ x) < m x [12], we have g lh (M) < g(M) < g uh {M), 

where 



1 / 2M 




even M 

(78) 

odd M 
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Fig. 14. Transmission capacity vs. path loss exponent of slotted ALOHA with OSTBC, and with r tr = 3 m, ft — 3 dB, 
e = 0.0001 and p = 1. 

From this, it is clear that, regardless of whether M is kept fixed or allowed to vary with N, the best achievable 
scaling is costbc( £ ) = © ^V«^. This result, combined with that for minimum -rate codes above, indicates that 
OSTBC can only achieve sub-linear scaling 10 , in contrast to the spatial multiplexing systems considered previously. 

We now consider the question of how to select M to maximize the transmission capacity scaling for large N. This 
is tantamount to maximizing the leading multiplicative factor g{M). For maximum-rate codes, we will focus on 
g up (M), since it serves as an accurate approximation for g(M). It is obvious that g np (M) is maximized by selecting 
M = 2 for M > 2. Thus, the key finding is that for OSTBC, either the Alamouti code or SIMO transmission should 
be employed. Although it is hard to obtain a simple design rule which specifies which of these two coding schemes 
will perform better under certain network conditions, we can gain insights by numerical analysis. In particular, 
Fig. 14 plots the transmission capacity achieved by OSTBC for different M. Results are presented for N = 3, 
and the codes used for M = 2 and M = 3 are the maximum-rate codes given by (15) and (36) respectively. The 
"Analytical" curves are based on (57), and are clearly seen to match with the "Numerical" curves, obtained by 
taking p large in the outage probability expression obtained by substituting the relevant parameters from Table I 
into (28), and subsequently solving for A(e) using numerical techniques. We see that the Alamouti code performs 
better than SIMO for a < 4.7, and vice-versa for a > 4.7. 

10 A similar scaling result was obtained in [12], based on the corresponding SINR approximation illustrated in Figs. 4-6. 
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C. Comparison 

For large N and assuming e > F su (/3) in (57), spatial multiplexing achieves a higher transmission capacity 
scaling than OSTBC, at high transmit SNR p, if 

W<#)(^ (1-#)) S for ZF 

[ +a> v 7 2 (79) 

^WX^ijtTfl-*)) 1 for MRC . 
Clearly, these conditions depend on the specific OSTBC employed. As expected, we see that by incorporating a 
higher rate code, whilst keeping all other parameters fixed, the possibility of the condition (79) being satisfied is 
decreased. Moreover, by noting that 7Vj(M) < M 3 for all practical OSTBC codes, the condition in (79) becomes 
more likely to be satisfied as the path loss exponent a increases; indicating that environments with higher path loss 
exponents are more beneficial for spatial multiplexing, compared with OSTBC; and vice-versa for low path loss 
environments. 

Figs. 15 and 16 show the transmission capacity achieved by spatial multiplexing with ZF receivers and OSTBC for 
different antenna configurations. The curves for both spatial multiplexing and OSTBC are based on (57). The results 
show that the relative transmission capacity of the transmission schemes depends on both the antenna configuration, 
the path loss exponent a and the SINR operating value. In particular, we observe that OSTBC performs best at small 
a (e.g., a = 2, corresponding to free space environments), whereas spatial multiplexing performs the best at large 
a (e.g., a = 6, corresponding to indoor environments). These observations agree with our analytical conclusions 
put forth above, based on (79). 

We also observe that ZF receivers are preferable over MRC receivers for a sufficiently high SINR operating value, 
i.e., p = dB, and vice-versa for a low SINR operating value, i.e., (3 = — 5 dB. This is because for high SINR 
operating values, the outage probability when there is no multi-node interference, i.e., the outage probability of 
the single user MIMO system, is higher for MRC than ZF, and offsets any positive gains resulting from the signal 
power that MRC has over ZF receivers. However, for sufficiently low SINR operating values, and as indicated by 
previous analysis, the outage probability when there is no multi-node interference is negligible for both receivers, 
and thus MRC is preferable over ZF. 

VI. Network Throughput Comparison between Slotted ALOHA and 

Coordinated Access protocol 

In this section, we compare the multi-antenna transmission schemes employing a slotted ALOHA MAC protocol, 
described in Section II, with a baseline single antenna transmission scheme employing a CA MAC protocol. It is 
worth noting that the CA protocol is an idealized protocol, since the overhead involved in achieving full coordination 
is prohibitive in practice for ad hoc networks. We will show that the simple multi-antenna slotted ALOHA MAC 
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Fig. 15. Transmission capacity vs. path loss exponent of slotted ALOHA with spatial multiplexing and OSTBC, and with 

N = 5, r tr = 1 m, f} = dB, N = 5, p = 30 dB, e = 0.15 and p = 1. 
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Fig. 16. Transmission capacity vs. path loss exponent of slotted ALOHA with spatial multiplexing and OSTBC, and with 

N = 5, r tr = 1 m, j3 = -5 dB, iV = 5, p = 30 dB, e = 0.15 and p = 1. 

protocol can in fact yield even higher throughputs than the tightly scheduled CA MAC protocol in various practical 
scenarios. This shows that not only does the use of multiple antennas allow for a decentralized MAC, but it can 
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actually also improve throughput in certain scenarios. We start by describing the CA protocol model. 

A. Coordinated Access Protocol Model 

For the CA protocol, transmissions are scheduled such that interference is minimized at each receiver. Since the 
strongest interferers are those closest to the receivers, we consider a CA scheduling protocol which removes the 
closest interferers within a guard zone around each receiver. Note that, although removing these strong interferers 
leads to a gain in throughput, this gain comes at the cost of significant overhead and synchronization requirements; 
something which is clearly undesirable in an ad hoc network setting. 

Scheduling nodes to minimize interference has been considered for a long time, however there have been few 
results where the transmitting nodes form a PPP. In [28], a circular guard zone of radius r gz around the typical 
receiver was considered, but neglected a guard zone around every other receiver in the network. Recently in [7], 
transmissions were scheduled only if the receiver has no interferers within a circular guard zone also of radius 
r gz . However, the throughput suffers when the network becomes dense, as the probability of a receiver with no 
interferers within its guard zone decreases. 

In this section, we consider an alternative model where transmissions are scheduled in a TDMA manner. As with 
other common CA protocols, the throughput of the TDMA scheme we consider does not go to zero as the network 
becomes dense. For analytical tractability, we assume each receiver is located at the center of a non-overlapping 
square of side length 2r gz , arranged according to a lattice structure. The transmitting nodes, distributed according 
to a PPP, which lie in the same square are each assigned a time slot of the same duration as the slots used in the 
slotted ALOHA protocol. We assume a dynamic slot allocation where the number of slots in a particular square is 
equal to the number of transmitters in that square. We note that although this may not be completely practical due 
to the deterministic placement of each receiver, this model correctly captures the fact that each receiver is spatially 
separated such that interference is minimized, and is very useful for comparison purposes. 

The CA protocol is illustrated in Fig. 17, where we show a network before and after scheduling. The scheduled 
transmitter transmits to the receiver at the center of the corresponding square during the assigned time slot. In Fig. 
17, the receivers are shown as crosses, while transmitters as dots. As shown, after applying the CA protocol, there 
is at most one transmitter in each square, which transmits to the receiver at the center of the square. 

Under this model, the following theorem presents a new closed-form upper bound for the network throughput. 

Theorem 3: The throughput of the CA MAC protocol employing single antennas is upper bounded by 



a \ yi j \a a \v 7rr tr 

(80) 

where Aca = (1 — e ~ 4Ar s z )/(4r gz ), and 2-F1O, ■; ■; ■) is the Gauss hypergeometric function [24, Eq. (9.6.2)]. 
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Fig. 17. Coordinated access protocol model: before and after scheduling with A = 3 m 2 . 

Proof: See the Appendix E. ■ 
Note that the same throughput can also be obtained by considering an alternative but similar system where the 
transmitting nodes around the square guard zone transmit with probability (1 — e^ 4Xr ^)/(4Xr^ z ) at each time slot. 
Thus the transmission probability of the CA protocol adaptively adjusts to the intensity of transmitting nodes. This 
is in contrast to the fixed transmission probability employed in slotted ALOHA. 

For the case of dense networks, the throughput of the CA protocol approaches a constant as A — > oo, and is given 
by substituting Aqa = AJ2? A = into (80). This constant value reflects the fact that the CA protocol ensures 

that only one node within a square of length 2r gz is allowed to transmit. This models, for example, a system using 
a TDMA protocol when a large number of transmitting nodes causes the throughput to reach a saturation point, 
and the only transmitting nodes are the ones allocated to the time slot corresponding to the current transmission 
period. 

Fig. 18 presents the throughput vs. the guard zone parameter r gz for different intensities. We see that the analytical 
upper bounds in (80) accurately match the Monte Carlo simulated curves in most cases. We also notice that there 
is an optimal value of r gz which maximizes the throughput of the system. In addition, we see that the dense 
network throughput, obtained by substituting \(? A into (80), closely matches the Monte Carlo simulated curves for 
transmitting node intensities as small as 0.03. 

B. Comparison 

In this section, we compare the slotted ALOHA MAC protocol employing multiple antennas with the CA MAC 
protocol employing single antennas. When considering the CA protocol we choose the guard zone parameter r gz 
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Fig. 18. Throughput of the CA protocol vs. guard zone parameter with p — 10 dB, r tr = 1.5 m, a = 4 and /3 = 5 dB. 

which maximizes the throughput for each value of A. Also, when considering the slotted ALOHA protocol we only 
utilize one transmit antenna and multiple receive antennas. Note that our analysis in Sections III and IV shows that 
utilizing more transmit antennas for spatial multiplexing and OSTBC may lead to a higher throughput in different 
networking scenarios. However, as our purpose is to show the practicality of the slotted ALOHA protocol, we 
consider the simple case when only one transmit antenna is utilized. Note that in this scenario, the performance 
of spatial multiplexing and OSTBC are the same. We assume both protocols use the same initial intensity of 
transmitting nodes A, thereby allowing for a fair comparison. 

Fig. 19 presents the intensities and transmission probabilities for which the throughput of slotted ALOHA is 
greater than for the CA protocol. Results are presented for different antenna configurations, and the respective 
throughput values are calculated from (24) and (80). Note that the throughput for slotted ALOHA can also be 
obtained from (24). We clearly see that increasing the number of receive antennas leads to an increase in the 
intensity range for which the slotted ALOHA throughput is greater than the CA throughput. 

In addition, Fig. 19 shows that slotted ALOHA outperforms the CA protocol for nearly all transmission probability 
values when the network is sparse. However, as the intensity increases, the transmission probability has to decrease 
in order for the throughput of slotted ALOHA to be greater than the CA protocol. Although not shown in Fig. 19, 
the throughput is an increasing function of the number of antennas. 
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Fig. 19. Comparison between slotted ALOHA and the CA protocol. The shaded area between the solid lines shows the region 
where the throughput of slotted ALOHA with M = 1 in (24) is greater than the CA protocol in (80); for N = 2. The dashed 
lines shows that the region expands for N = 3, and the dotted lines shows further expansion for iV = 4. Results are shown for 
M = 1, p = 10 dB, TV = 2 m, a = 3 and j3 = 5 dB. 

VII. Conclusion 

Open-loop point-to-point transmission schemes and slotted ALOHA MAC are practical choices for ad hoc 
networks, due to their relatively low feedback requirements and decentralized structure, respectively. Within this 
framework we have analyzed important network performance measures for spatial multiplexing transmission with 
MRC and ZF receivers, and for OSTBC, using tools from stochastic geometry. We derived new expressions for the 
outage probability, which were subsequently used to evaluate the network throughput and transmission capacity. 

Based on our analysis, we demonstrated that the optimal number of transmit antennas, and the relative performance 
of the spatial multiplexing and OSTBC systems is dependent on a number of system parameters, including the chosen 
SINR operating value and node intensity. In particular, from a throughput perspective, our analysis has led to the 
following design guidelines: 

• For dense networks and high SINR operating values, cyclic antenna diversity codes utilizing the maximum 
number of transmit antennas is the preferred scheme. 

• For low SINR operating values, spatial multiplexing with MRC receivers, while transmitting with the maximum 
number of data streams, is the preferred scheme. 

For applications operating under a strict outage constraint, the transmission capacity is a preferable performance 
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measure. From a transmission capacity perspective, the following were observed: 

• For sufficiently high path loss exponents, spatial multiplexing is preferred over OSTBC, and vice-versa for 
low path loss exponents. 

• In the large antenna N regime, and if the self-interference is sufficiently smaller than the SINR operating 
value: 

- MRC and ZF receivers can achieve linear scaling, while OSTBC can only achieve sub-linear scaling. 

- For OSTBC, either the Alamouti scheme or SIMO is the optimal scheme, depending on the path loss 
exponent and number of receive antennas. 

Finally, we demonstrated the practicality of the multi-antenna slotted ALOHA MAC system in ad hoc networks, by 
comparing it with a baseline single antenna CA MAC system. We showed the interesting result that in many practical 
scenarios, the multi-antenna slotted ALOHA MAC system can actually outperform the idealized single antenna CA 
MAC in terms of network throughput. Further, we demonstrated that increasing the number of antennas can have 
the benefit of increasing the range of node intensities for which the use of slotted ALOHA MAC outperforms the 
idealized single antenna CA MAC. 
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Appendix 

A. Derivation of SINR Approximation for OSTBC 
We require the following lemmas: 

Lemma 1: Let h ~ CNnxi (0jvx1,Ijv) and g ~ CNnxi (OjvxIjIjv) be independent random variables. Then 

~Gamma(0,l) (81) 

is independent of h. 

Proof: See [29]. ■ 
Lemma 2: Let h ~ CNnxi (Ojvxi, Ijv)> 
g <~ CMnxi (0jvx1)Iat)> x <~ CAf(0, 1) and y <~ CAf(0, 1), all mutually independent. Then 

|htg| 2 < st Ihtg + xyl 2 (82) 

where 

A < st B => F A (z)>F B (z), forzeR+, (83) 
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for random variables A and B with c.d.f.s F^(-) and Fb(-) respectively. 

Proof: Setting X\ = h^g, the distribution of X\ conditioned on h is CM (0, ||h|||-). Thus, \Xi\ 2 conditioned 

on h has the same distribution as ||h|||.Yi, where Y\ ~ Exp(l). Similarly, setting X 2 — h^g+rry, \X 2 \ 2 conditioned 

on h has the same distribution as (||h||f, + \x\ 2 ) Y 2 , where Y 2 ~ Exp(l). (82) now follows. ■ 
Lemma 3: Let Xi ~ Gamma (kg, 6), for i = 1, . . . , L, all mutually independent. Then 

X t ~ Gamma k ^ 6 ) ■ ( 84 ) 

i=l \i=l / 

Proo/- See [30]. ■ 
Lemma 4: Let X x ~ Gamma(fci,0) and X 2 <~ Gamma(fc 2 , 6*). If fci > k 2 , then 

E [Xf] > E [X%] (85) 

for any n € Z + . 

Proof: The nth moment of X ~ Gamma (k, 9) is 

E[x] - r(fc) • (86) 

The proof follows upon showing the positivity of the derivative of (86) w.r.t. k. ■ 

Lemma 5: Let h ~ CAf Lxl (0 Lxl ,I L ), g ~ CAf Lxl (0 Lxl ,I L ) and y ~ CAf Tx i (0 Txl ,I T ), all mutually 
independent. Then Z — ||h|p+||y|| 2 nas approximately the same distribution as 

Z~ ~ Gamma f -, 1 ) . (87) 



\L + T 

Proof: We justify this approximation by first noting that E[Z] = E[Z~]. Also, applying Lemma 1 and Lemma 
4 to Z yields: Z lb < Z < Z uh , where Z lh ~ lim ? ^ Gamma (<r, 1) and Z ub ~ Gamma (1, 1). Now, Z Z lb 
when L -> and Z -> Z ub when L^L + T^T^O, implying that Z -> Z~ at the two extremes: -> 
and j-q^ — > 1. For values away from these two extremes, our approximation is motivated by noting that for a fixed 
L + T, applying Lemma 2 to Z and Lemma 4 to Z~ reveals, respectively, that Z and Z~ both increase with L. ■ 
We are now in a position to derive our approximation for the OSTBC SINR. We first outline the main issues 
which make dealing with the exact SINR (19) difficult, and provide some justification for our approximations which 
are used to overcome these. 

• The first challenge is caused by mutual dependencies between random variables in the SINR. In particular, 
the normalized interference power, tCe,$2> is comprised of a sum of dependent random variables. As in [12, 
23], we will neglect this mutual dependence, which is justified in [23] by noticing that the common terms in 
JCe t j2 are multiplied by independent random variables, so they should be "nearly independent". In addition, the 
numerator and denominator in (19) are dependent. Once again, as in [23], we neglect this dependence, with 
the argument that there are sufficient independent terms in the numerator and denominator that the dependence 
is weak. 
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• The second challenge is the difficulty in computing the exact marginal distribution of /C^y- analytically, even 



under the independence assumption discussed above. This appears intractable, and to address this we seek an 
accurate approximation by fitting a gamma distribution. Such distributions are known to be exact for various 
cases and are expected to be quite accurate in general. They are also analytically friendly, allowing us to 
invoke the useful properties given in the lemmas above. 
To most clearly illustrate our methods, we will first present our derivations for a specific OSTBC code; namely, 

the OSTBC code in (35). We will then give a brief discussion to highlight how the derivation extends to general 

OSTBC codes. It is convenient to define the operator vecdim(-), which takes as input the inner product of two 

vectors with equal dimension, and returns the dimension of these vectors. 

1) Derivation for the OSTBC Code in (35): Without loss of generality, we focus on decoding x ,i- The 

£th row of ?fe(V) and are given respectively by 





~Vi,i\ and 



(88) 



Substituting (35) and (88) into (18), the data estimate for x 0j i can be written as 




(89) 



where 



(90) 



with 




(91) 



N 




fe=l 



N 




(92) 



fe=i 



N 




N 




(93) 
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Note that the We, q s and Zi, q s can each be interpreted as the inner product of two equal-length vectors; in each 
case, with one vector having elements drawn from either ±Ho or ±Hq, and the other vector having elements drawn 
from either ±H f or ±Hj. For example, We,2 in (92) can be written as 

We, 2 = w • wg (94) 
where w and are 2N x 1 vectors given respectively by 

w = [^0,1,1; _ ^0,1,2, ••■ , ^0,jV,l> — /io,W,2] T (95) 

and 

= [he,i,2, • • • , he,N,2, h* t ,N,i\ T ■ (96) 

With this interpretation, it is clear that vecdim(T4^i) = AN and vecdirr^W^) = vecdim(^ j2 ) = vecdim(Vl^ ! 3) = 
vccdim(Z £:3 ) = 2N. 

From (89), the SINR can be written as 

p IITJ 112 

RMr° 1 1 1 1 F 

70STBC.1 = p ^ 1 p — 7 (97) 

RM l^D e £<S> \D e \« K -i,'E + 1 

where K, ey - = E ''- 1, ',?£ ? v f [hI ' ] ■ Note that 

t,z - ll H o|| F 

ICe,j2 — + ^,2 + fce,3 (98) 



where 



\We, 



K <<' " w (99) 



and, for q = 2, 3, 



E^ 



I 1 2 

\We, q xe, q + Ze tq xe iq \ 
IIHol 12 



if 

|W^| 2 + |Z,, 9 | 2 



(100) 



ll n o|| F 

It is clear that all three terms in (98) are mutually dependent. Moreover, there is dependence within each term, 
since \Wg, q \ 2 , \Zi tq \ 2 , and ||Ho|||. involve common elements. As stated previously, we neglect these dependencies. 

I W I ^ 1 7 I ^ 

Thus applying Lemma 5 to pj^yp and , followed by applying Lemma 3 reveals that K,e, q is approximately 

distributed as 

Gamma (vecdim(W^ g ) + vecdim(i^ ;9 ), 1) (101) 
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for all q = 1,2, 3. (Note that for the particular code considered here, Z^i — 0.) Now, from (35) and (90), it can be 
seen that 

vecdim(W e , g ) + vecdim(Z e , q ) = NNj , (102) 

5=1 

and we can therefore deduce that /C^y; is approximately distributed as Gamma (Ni/M, Ni). Finally, with the 
aforementioned assumption that the numerator and denominator of the SINR (97) are independent, we obtain the 
desired approximation. 

2) Extension to General OSTBC Codes: Here we give some details of how to generalize the previous 
derivation for arbitrary OSTBC codes. Since this extension is quite straightforward, the discussion is kept brief. 
The challenge, once again, is to approximate the distribution of fCe s = Ex *[^*®^<! x i)lli] m ^ & before, it 

II^oIIf 

can be shown that 



^£.s — y^^e, q , K-i.q — - — £ ' < ?J rj .I e ''^ (103) 

where W^_ q and Z^ q are each composed of the inner product of two equal-length vectors; one vector having elements 
drawn from ±Ho or ±Hq without repetition, the other having elements drawn from ±H^ or ±Hj without repetition. 
As evident from the previous example, the specific set of selected elements as well as the dimension of the vectors 
comprising the inner products are dependent on the particular OSTBC code employed. 

Under the same independence assumptions as before, we apply Lemma 5 to , and ,Lf ' ", 2 , followed by 

II^oIIf II^oIIf 

applying Lemma 3, to find that JCi t<l is approximately distributed as 

Gamma (vecdim(I^ j9 ) + vecdim(Z^ 9 ), 1) . (104) 

It follows that /Cf,s is approximately distributed as Gamma (Nj/M, 1). The desired SINR approximation is then 
obtained by invoking the assumption that the numerator and denominator and independent, as before. 

B. Proof of Theorem 1 

We first rewrite the c.d.f. of 7 as 



W 

■m<j-Y-l 



A Pr(A) 



= 1 - Pr ( \Xi\- a ^m < y - Y - 1 A j Pr(A) - Pr ( ]T \X t \- a * t 
= 1 - Pr ( ]T \Xe\- a V i0 <^-7-lA Pv(A) (105) 



where A denotes the event ^ — Y > 1 and A denotes the complement of A. Applying [6, Eq. (113)] along with 
some algebraic manipulation, the c.d.f. of 7, conditioned on W and Y, can be written as 

F M yW) = 1 - E V ,, r , ^ ' ((W -Y0-(i)-E [*=] ) Pr(A) 



k\r(i-m 

fe=0 V a ) 
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where ^ ~ "J^o f° r all Averaging with respect to (w.r.t.) W gives 



~ (-7rpAT(l-|)^E[*H]) 



fc=0 



fc!r(i - f ) 



(106) 



where 



POO 

li= / (w-Yfi-P)-* f w (w)dw 

g(y+jj m-l 



e » 
r(m)^ 



E 



fcO 



m — 1 



2fc 



(/?(y + i))T [ m -£- ]ff 



771 — t— ±zr 



(107) 



Substituting (107) into (106), we have 



gg+i) m-l 



F 7 ,y(/3) - 1 - 



E 



m-l\ ^/3(y + l) 



< oo 

E 



/ _ n^(i-g)^E[,a ] y r(m _ / _^ ) 



(108) 



r(m) t J \ o J fcir (i - 2*) 

The remaining challenge is to remove the infinite series in (108). To this end, we apply the identity [24, Eq. (24.1.3)] 

r(m-^-f) 



r(i-f) 

and Dobihski's Formula [31], which gives 

7rpAr(l-i)/3S E^* = J 



m-t-i /OlA ' 

= ( _ ir -,-l ^ s(m _^ +1) W 

i=0 ^ ' 



(109) 



m-l 



r(m) 



E 



m-t-1 



J2 s(m-e,i + i)l- \ E^i) 

i=o j=0 



m-l\ / /?(y + i) 



TrpAT (l - £) 0£e 



(HO) 



Averaging w.r.t. F, we have 



F 7 (/3) 



= 1 (-ir-^ 



r(m) 



E 



m — 1\ / /? 



^ s(m-£,i + l) (111) 



i=0 



E 5 (^') 

j=0 



7rpAT(l-£)^E[*i 



where 



/•CO 

X 2 = / e-^ (y + lfi Y (y)dy 
Jo 



1 

= r(u)T u S 



T=0 



Finally, by noting that 



E 



^7 



A 1 
+ T 



^°r(n + f) 

I» 



T(t + u) . 



(112) 



(113) 



and substituting (112) and (113) into (111), we obtain the desired result. 
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C. Proof of Corollary 5 

We first note that by taking a Taylor expansion of (28) around A = 0, and then finding the inverse of the resulting 
expression w.r.t. Xp, the transmission capacity can be written in the general form 



c(e) = 

ir(3irlE g (SNR ) 



where 



ff (SNR )= ^ {w-YP-PY" fw(w)dw . (115) 

J0(Y+1) 

To proceed, we note that if X <~ Gamma(n, then lim^oo X ^!f ~ CA/ r (0, 1). Now if lim^oo y 7 ™? — > and 
lim„_ ) . (X) n$ — > c where c is a constant, then lim^oo X — > c. We see from Table I that as N — > oo with M = kJV 
where < k < 1, -> 0, and -> 0, and thus 

W^mO, F^uT and * -> nft . (116) 

The result follows by substituting (116) into (114). 

D. Proof of Proposition 1 

Let M = f(N). For large M and iV, 



/(AO 



where 



l + f orZ F 

< = < " • (118) 

forMRC 

Since M < N, f(N) must satisfy f(N) < N. Therefore, f(N) must scale linearly or sub-linearly with N. If 
/(AT) is linear in N, then we know that the transmission capacity is linear in N also. However, if f(N) scales 
sub-linearly, then the transmission capacity (117) becomes 

C = e(iv^/ 1 -^(7V)) (ii9) 

which is also sub-linear in N. 
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E. Proof of Lemma 3 

The intensity Aca after applying the CA protocol is obtained by noting that at most one node is allowed to 
transmit in each square of length 2r gz , and 1 — e b z is the percentage of receivers with at least one node inside 
the square guard zone. By denoting the new aggregate interference after applying the CA protocol as pX$ CA , and 
noting that the typical transmitter-receiver channel, \ho\ 2 , is exponentially distributed, we can write the throughput 

as 



where ( — and /z #CA (x) is the p.d.f. of <1>ca- Note that the third line in (120) was obtained by applying 
standard results in stochastic geometry (see e.g., [16, pp. 231]). Note that the integral in (120) is over the entire 
infinite plane, except for a square guard zone centered at the origin of length 2r gz . Due to the difficulty in evaluating 
an exact closed-form solution to this integral, we obtain an upper bound by assuming a circular guard zone of radius 
2 r gz/Vn centered at the origin. The use of a circular guard zone as an upper bound can be seen by first noting 
that the same number of interferers, on average, are removed from within both the square and circular guard zones. 
Second, by denoting 5* as the interferers which do not fall within both guard zones, then the interferers in S which 
are enclosed by the circle of radius 2r gz /y / 7r are closer to the origin than the interferers in S which are enclosed 
by the square of side length 2r gz . This implies that the use of the circular guard zone removes interferers which are 
closer to the typical receiver, than the interferers which are removed using the square guard zone. The throughput 
upper bound is thus given by 



Finally, we solve the integral in (121) using [32, Eq. (3.259.2)-(3.259.3)] followed by some algebraic manipulation. 
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